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Abstract 



^ ■ The prospects for the measurement of the pseudoscalar admixture in the Htt 



coupUng to a Standard Model Higgs boson of 120 GeV mass are discussed in a 
quantitative manner for e~^e~ collisions of 350 GeV center-of-mass energy. Specific 
angular distributions in the h r+r^, p^i?^[i>r) decay chain can be used to 

probe mixing angles of scalar-pseudoscalar Htt couplings. In the discussion of the 
feasibility of the method, assumptions on the properties of a future detector for an 
e'^e~ linear collider such as TESLA are used. The Standard Model Higgsstrahlung 
production process is taken as an example. For the expected performance of a typical 
Linear Collider set-up, the sensitivity of a measurement of the scalar-pseudoscalar 
mixing angle turned out to be 6°. It will be straightforward to apply our results 
to estimate the sensitivity of a measurement, in cases another scenario of the Higgs 
boson sector (Standard Model or not) is chosen by nature. The experimental error 
of the method is expected to be limited by the statistics. 
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1 Introduction 



The transverse spin effects in r pair production can be helpful to distinguish between the 
scalar J^*^ = 0^+ and pseudoscalar J'^'^ = 0^+ natures of the spin zero (Higgs) particle, 
once it is discovered in future accelerator experiments. To address resolution issues, it 
is necessary to perform Monte Carlo studies, where the significant details of theoretical 
effects and detector conditions can be included. To enable such studies we have extended 
the algorithm of Refs. [1,2] of the TAUOLA r-lepton decay library [3-5] to include the 
complete spin effects of r leptons originating from the spin zero particle. 

In Refs. [6,7] the reaction chain e+e" ^ Z{H/A°), H/A^ t+t', r± n^M^r) 
was studied. It was found that even small effects of smearing seriously deteriorate the 
measurement resolution. However, using the TAUOLA spin interface, we devised a very 
promising method for the measurement of the Higgs boson parity, see Ref. [8]. It turns out 
that the spin effects of the decay chain H/A^ t~^t~ p~^u^p~i'^ Tr+vr^z/^TT^Tr^z/^ give 
a parity test independent of both model (e.g. SM, MSSM) and Higgs boson production 
mechanism (e.g. Higgsstrahlung, WW fusion). In the rest frame of the p'^p~ system 
we defined the acoplanarity angle (p* as the one between the two planes spanned by the 
immediate decay products (the tt^ and 7r°) of the two p's. This angular distribution of the 
T decay products, which is sensitive to the Higgs boson parity, once additional selection 
cuts are applied, is measurable using typical properties of a future detector at an e^e~ 
linear collider. Using reasonable assumptions about the SM production cross section and 
about the measurement resolutions we have found that, with 500 fb~^ of luminosity at a 
500 GeV e+e" hnear collider, the CV of a 120 GeV Higgs boson can be measured to a 
confidence level greater than 95%. 

In Ref. [9] we demonstrated that a measurement of the r impact parameter in one- 
prong T decay is useful for the determination of the Higgs boson parity in the H/AP — > 
T+T~; — > p^Dt{vt-) decay chain. We estimated that, for a detection set-up such as 
TESLA, use of the information from the r impact parameter can improve the significance 
of the measurement of the parity of a Standard Model 120 GeV Higgs boson to ~ 4.5(T 
and in general by a factor of about 1.5 with respect to the method where this information 
is not used. So far we have not exploited the possibility of using decay modes other than 

— > p'^i/t{vt). Additional modes are expected to further increase the separation power. 

In this paper we study the more general case where mixed scalar and pseudoscalar 
couplings of the Higgs boson to r leptons are simultaneously allowed, see e.g. Ref. [10]. 

Our paper is organized as follows: In Section 2 we present basic properties of the 
density matrix for the pair of r leptons produced in Higgs boson decay. In Section 3 we 
define our observable and in Section 4 our Monte Carlo set-up. Our results are presented 
in Sections 5 and 6, first with an idealized detector set-up and then with more realistic 
assumptions on the detector and integrated luminosity. A summary. Section 7, closes the 
paper. 
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2 Spin weight for the mixed scalar— pseudoscalar case 



Let us here, only very briefly describe the basic properties of the spin correlations and 
their implementation in our Monte Carlo algorithm. We will not repeat the detailed 
description of the method (which can be found in Ref. [3]) or the algorithm (which is 
given in Ref. [6]). We will discuss the points necessary to understand the case of mixed 
scalar-pseudoscalar coupling of hrr. 

The main spin weight of our algorithm for generating the physical process of r lepton 
pair production in Higgs boson decay, with subsequent decay of r leptons as well, is given 
by 

^^ = i(l + EE^^^-^i^2), (1) 
V i=i j=i J 

where hi and h2 are the polarimeter vectors that depend respectively on decay products 
momenta; Rij is the spin density matrix. For the mixed scalar-pseudoscalar case, when 
the general Higgs boson Yukawa coupling to the r lepton 

f (a + z675)r (2) 

is assumed, we get the following non-zero components of Rij: 

R - 1 n -n -'ll^ n - H _ .o^ 
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where /3 = yl — If we express Eq. (jSJ with the help of the scalar-pseudoscalar 

mixing angle 0: 

f A^(cos + i sin 075)r, (4) 
the components of the spin density matrix can be expressed in the following way: 
_ _ cos 02 - sin 02 D _ D _ 2cos0sin0/5 

-Kll — -^22 — ,0 /32 I '■ 12' -"12 — --K2I — ,0 02 , '■ 72' 

COS 0^ + sm 0^ COS 0^ + sm 0^ 

In the limit /3 — > 1 these expressions reduce to the components of the rotation matrix for 
the rotation around the z axis by an angle —20: 

R\\ = -R22 = cos 20, i?i2 = — -R21 = sin 20. (6) 

The Higgs boson parity information must be extracted from the correlations between 
r"*" and r~ spin components, which are further reflected in correlations between the r 
decay products in the plane transverse to the r+r" axes. The same will now apply to 
the mixing scalar-pseudoscalar case. To better visualize the effect to be measured, let us 
write the decay probability for the mixed scalar-pseudoscalar case, using the conventions 
of Ref. [11]: 

r(/i™. ^ r+r-) ~ 1 - s\\s\ + s^i?(20) s\ , (7) 

where -R(20) can be understood as an operator for the rotation by an angle 20 around the 
II direction. The and are the polarization vectors, which are defined in their 
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respective rest frames. The spin quantization axes are oriented in the r~ flight direction. 
The symbols ||/_L denote components parallel/transverse to the Higgs boson momentum 
as seen from the respective rest frames. 

It is straightforward to see that the pure scalar case is reproduced for (j) = 0. Then 
Rii = +1, R22 = +1 and R^s = — 1 are obtained, and the limit P ^ I does not need 
to be taken. For = 7r/2 we reproduce the pure pseudoscalar case. We get -Rn = — 1, 
R22 = — 1 and i?33 = — 1. Also in this case, the /3 — > 1 limit was not needed. 

3 The acoplanarity of the and p~ decay products 

To facilitate reading, let us recall here some elements of the observables that were pre- 
sented in Refs. [8, 9] and can be used to measure the Higgs boson parity. We will stress 
only those points that required modification. The method relies on measuring the acopla- 
narity angle of the two planes, spanned on p"^ decay products and defined in the p'^p~ pair 
rest frame. For that purpose the four-momenta of vr^ and vr'' need to be reconstructed 
and, combined, they will yield the p^ four-momenta. All reconstructed four-momenta 
are then boosted into the p'^p~ pair rest frame. The acoplanarity angle ip*, between the 
planes of the p"*" and p~ decay products is defined in this frame. In the previous papers 
only the range vr was interesting and thus reconstructed, as this was sufficient to 

distinguish between two possibilities: scalar or pseudoscalar Higgs boson, differing by the 
sign of the transverse spin correlation. The angle was defined with the help of its cosine 
and with the help of the two vectors n± normal to the planes namely n± = p^± x p^o, 

COS 09 = ^=11 r. 

r- |n+||n_| 

For the present use, such a definition is insufficient. As can be seen from Eq. (0) the 
correlation, in the case of the Higgs boson of combined scalar and pseudoscalar couplings 
of Eq. (0} and the mixing angle (p, is between transverse components of r+ spin polar- 
ization vector and transverse components of t~ polarization vector rotated by an angle 
2(j). Therefore the full range of the variable 27r is of physical interest. To 

distinguish between the two cases (p* and 27r — ip* it is sufficient, for example, to find the 
sign of Ptj- ■ n+. When it is negative, the angle p* as defined above (and in the range 
< p* < 7i) is used. Otherwise it is replaced by 2it ~ p>*. If no separation was made, the 
parity effect, in case of mixed hrr coupling, would wash itself out (see Fig. |2l later in the 
text). For the graphical representation of the definition of the angle p*, see Fig. ^ The 
figure visualizes the relation between the observable and Eq. ((7j) as well. 

Additional selection cuts need to be applied. Otherwise the acoplanarity distribution 
is not sensitive to transverse spin effects (and thus to Higgs boson parity) at all. The 
events need to be divided into two classes, depending on the sign of yiy2, where 

Ej^+ + E^o Et^" + E^o 

The energies of 71^,71^ are to be taken in the respective rest frames. In Refs. [8,9] 
the methods of reconstruction of the replacement rest frames were proposed with and 
without the help of the r impact parameter. We will use these methods here as well, 
without any modification. 
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Figure 1: Definition of the p'^p~ decay products' acoplanarity distribution angle (p* , in 
the rest frame of the p^ p" pair. The range for ip* is < if* < 27i. Note that, for better 
visualization, we use in this figure the momenta of ir^ and p^ (rather than n^'s from p^ 
decays) to define the planes. The two ways of defining the planes are equivalent if no 
reconstruction errors are taken. 

4 The Monte Carlo 

If any non-zero CP-odd admixture to the Higgs is present, not only is the distribution of 
the Higgs decay products modified, but also the distribution of its production angle [10- 
12]. In this paper, we simulate production angular distributions as in the SM, but this 
assumption has no influence on the validity of the analysis. In particular, the detection 
efficiencies for pure CP-even and pure CP-odd Higgs bosons do not differ significantly. 
In order to study the sensitivity oi h t^t^ observables, we assume a production rate 
independent of the size of the CP-odd admixture, i.e. the SM production rate of a CP-even 
Higgs. 

The production process e+e" — >■ Zh — > p'^p~{q(l)T~^T~ has been chosen, as an repre- 
sentative example, and simulated with the Monte Carlo program PYTHIA 6.1 [13]. The 
Higgs boson mass of 120 GeV and a center-of-mass energy of 350 GeV was chosen. The 
effects of initial state bremsstrahlung were included. For the sake of our discussion and 
in all of our samples the r decays have been generated with the TAUOLA Monte Carlo li- 
brary [3-5]. As usual, to facilitate the interpretation of the results, bremsstrahlung effects 
in decays were not taken into account. Anyway, with the help of additional simulation, we 
have found this effect to be rather small. To include the full spin effects in the h — > t~^t~ , 
t"^ —>■ p^^'t-(z/t-), Tc'^TT^ decay chain, the interface explained in Ref. [6] was used, with 

the extensions discussed in Section 2. 
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5 Idealized results 



5.1 Resolution parameters 

To test the feasibihty of the measurement, some assumptions about the detector effects 
had to be made. We include, as the most critical for our discussion, effects due to 
inaccuracies in the measurement of the tt^,7i^ momenta and of the leptons impact 
parameters. We assumed Gaussian spreads of the measured quantities with respect to 
the generated ones and we used the following algorithm to reconstruct the energies of 
tt's in their respective rest frames, exactly as in the case of the studies presented in 
Refs. [8,9]. 

1. Charged-pion momentum: We assume a 0.1% spread on its energy and direction. 

2. Neutral-pion momentum: We assume an energy spread of . For the 9 

and angular spread we assume I^q- These 7r° resolutions can be achieved with 
a 15% energy error and a 27r/1800 direction error in the gammas resulting from the 
7r° decays. These resolutions have been approximately verified with SIMDET [14], a 
parametric Monte Carlo program for TESLA detector [15], as well as with other 
studies, see e.g. Refs. [16,17]. 

3. The reconstructed Higgs boson rest frame: We assume a spread of 2 GeV with 
respect to the transverse momentum of the reconstructed Higgs boson momentum, 
and 5 GeV for the longitudinal component, to mimic the beamstrahlung effect. 

4. The impact parameter: The angular resolution of the r impact parameter 
has been simulated for a TESLA-like detector. The simulation is based on the 
anticipated performance of a 5-layer CCD vertex detector, as described in Ref. [15]. 
For Higgsstrahlung events with hsM t^t^ and — > p='=z7^(z/^) at m^gj^j = 
120 GeV and = 350 GeV, the angular resolution has been found [9] to be 
approximately 25°. 

5.2 Numerical results 

We have used the scalar-pseudoscalar mixing angle = f and, as the reference, we have 
used the pure scalar case = 0. 

In Fig. 121 the acoplanarity distribution angle (p* of the p^p^ decay products which 
was defined in the rest frame of the reconstructed p^p^ pair, is shown. Unobservable 
generator-level rest frames are used for the calculation of selection cuts. The two plots 
represent events selected by the differences of tt^tt^ energies, defined in their respective 

rest frames. In the left plot, it is required that yiy2 > 0, whereas in the right one, 
events with yiy2 < are taken. This figure quantifies the size of the parity effect in an 
idealized condition, which we will attempt to approach with realistic ones. The size of the 
effect was substantially diminished when a detector-like set-up was included for rest 
frames reconstruction as well, see Fig. 13 in exactly the same proportion as in Ref. [8]. 
The general shape of the distributions remained. 
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At the cost of introducing cuts, and thus reducing the number of accepted events, we 
could achieve some improvement of the method, as in Ref. [9]. If we require the signs 
of the reconstructed energy differences yi and 1/2 (Eq. (jS))) to be the same whether the 
method is used with or without the help of the r lepton impact parameter, only ~ 52% of 
events are accepted. The relative size of the parity effect increases. Results are presented 
in Fig. El 




0.25 ■ N„ r 




0.25 -N^ r 



Figure 2: The acoplanarity distribution (angle (f*) of the p~ decay products in the rest 
frame of the p~ pair. Gaussian smearing of it's momenta as described in Section 5.1 
is included. However, generator level rest frames are used. The thick line corresponds 
to a scalar Higgs boson, the thin line to a mixed one. The left figure contains events with 
?/i2/2 > 0, the right one is for yiy2 < 0. In our paper, that is for the 350 GeV e~*"e~ CMS 
(scalar 120 GeV mass) Higgsstrahlung production we took = 62.7-10^^ [fb] for the scale 

of the plot. In general case N„ = -^Ctotaii,^^ ^~ XH)BTZ{H — > t^t^){BTZ{t Pt^t)Y 
is a suitable choice. 



6 Simulation with detector effects 

In order to assess the possibilities for a measurement of the acoplanarity distribution 
described in Section 2, we perform a detailed simulation of Higgs bosons produced in the 
Higgsstrahlung process using PYTHIA 6 . 2 [13] for the production process and the modified 
version of TAUOLA described above to generate samples of signal events. These events 
are then passed through a simulation of the TESLA detector (SIMDET [14]) accounting 
for the acceptance and anticipated resolution of the tracking devices and calorimeters 
corresponding to the detector proposed in the TESLA TDR [15]. 

Signal samples ^ of 1 ab^^ at 350 GeV center-of-mass energy were generated for scalar- 

^Note that this integrated luminosity is larger by a factor of 2 than the one used in Refs. [8, 9] to 
estimate the sensitivity of our Higgs boson parity observable. Also, the Higgstrahlung production cross 
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Figure 3: The acoplanarity distribution (angle ip*) of the p~ decay products in the rest 
frame of the p'^p~ pair. Gaussian smearing ofn 's and Higgs boson momenta, as described 
in Section 5. 1 is included. Only events where the signs of the energy differences yi and y2 
are the same, if calculated using the method described in Ref. [8] and if calculated with the 
help of the r impact parameter Ref. [9], are taken. The thick line corresponds to a scalar 
Higgs boson, the thin line to a mixed one. The left figure contains events with yiy2 > 0, 
the right one is for yiy2 < 0. In our paper, that is for the 350 GeV e~^e~ CMS (scalar 
120 GeV mass) Higgsstrahlung production we took = 62.7 ■ 10~^ [fb] for the scale of 
the plot. In general case = -^(Jtotaii.^^^' XH)BTZ{H — >• t^t^)[BIZ{t pi^r)) is 
a suitable choice. 



pseudoscalar mixing angles = 0, vr/S and 7r/4. With detector simulation the r leptons 
decaying to tt^tt^ from Higgs decays were reconstructed as isolated jets with only one 
charged track (the reconstructed vr^) and additional neutral clusters (the reconstructed 
7r°). The vr^ and 7r° momenta were combined to form a reconstructed p^. The acopla- 
narity angle if* was calculated in the reconstructed p^p^ rest frame. Two event classes 
are formed according to the sign of yiy2, where yi and y2 are calculated in the labora- 
tory frame. The resulting ip* distributions for the three cases are shown in Fig. 0] as 
histograms, each containing about 0.5 ab~^ statistics. 

To extract the scalar-pseudoscalar mixing angles the functions a * cos{(f* — 20) + b 
(for yiy2 > 0) and a * cos(v?* — 20 + tt) + 6 (for yiy2 < 0) were used to fit 20 to the 
reconstructed acoplanarities y?*, gained from simulated detector signals. The constants a 
and b were additional free variables of the 3-parameter fit. The resulting functions are 
also shown as lines in Fig. 0] 

In order to assess the expected accuracy and a possible experimental bias of the 
measurement, the above procedure was repeated 400 times with acoplanarity distributions 

section (see e.g. [18]) is more than 2 larger at 350 GeV than at 500 GeV center-of-mass energy. On the 
other hand, here we do not use the information from the r impact parameter, which can be useful to 
improve the sensitivity of a measurement of the mixing angle 4>. 
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Figure 4: Distribution of the reconstructed acoplanarity angle (p* for </) = (full his- 
togram), (f) = 7r/8 (dashed histogram) and cj) = n/A (dotted histogram) for yiy2 > (left) 
and yiy2 < (right). The lines indicate the results of the corresponding fits (see text). 

extracted from independent samples of 1 ab~^ luminosity each, with a nominal value of 
(f) = 7c/4. Unlike what was done before in Fig. |3J the data for the two ranges of value of 
2/1^/2 were appropriately combined into one if* distribution before the fit. The new value 
of (f* for the case of yiy2 < had to be redefined as + vr for < < tt and ip* — n 
for IT < (f* < 271. The distribution of the fit results on 20 for each of the experiments is 
shown in Fig. El The mean value is 1.627 ± 0.014, compared to the | input value. The 
resulting bias of approximately 3° can probably be corrected in the future. The expected 
error on 20 is obtained as the width of this distribution. It amounts to 0.20 ± 0.01 rad, 
or approximately 12°. Thus, a precision on of approximately 6° can be anticipated for 
a SM Higgs cross section and h — >• r+r~ branching ratio at y/s = 350 GeV and 1 ab~^. 
Note that so far backgrounds neither from other Higgs boson decays nor from other SM 
processes have been considered. While previous studies [19] have shown that h t^t~ 
events can be selected without large backgrounds, some small deterioration and a further 
lowered signal efficiency are to be expected. Because of the small observed bias, it is not 
expected that systematic effects will limit the resolution even for production cross sections 
a few times larger than in the SM. 



7 Conclusions 

We have found that for an integrated luminosity of 1 ab~^, at 350 GeV center-of-mass 
energy, a high precision LC detector such as the proposed TESLA, should be able to 
measure the scalar-pseudoscalar mixing angle for the Htt coupling with 6° accuracy in 
the case of a Standard Model Higgs boson of 120 GeV mass. The experimental error is 
expected to be dominated by statistics. 
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However, if the production mechanism of the Higgs boson happened to be non-SM and 
larger, the systematic errors not studied so far and possible new, unknown phenomena 
may have some significant influence. On the basis of studies performed to date, we believe 
however that, if the cross section were somewhat larger than the Standard Model one and 
thus the uncertainty on the mixing angle due to statistics were not smaller than 4°, we 
do not expect the systematic error to be a problem. In the case of Higgs boson scenarios 
predicting even higher rates of observed h tt samples, the issue of the systematic 
error definitely needs to be re-addressed before any conclusion on measuring the scalar- 
pseudoscalar mixing angle in the hrr coupling with higher precision can be attempted. 

Finally, let us note that this method can be applied to measure the parity properties 
of other scalar particles, not necessarily only Higgs boson(s). 
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Figure 5: Distribution of the fit values for 20 from 400 independent samples, each corre- 
sponding to a luminosity of 1 ab~^, for a generated value of (p = ir/A. The curve represents 
the fit of a Gaussian to this distribution. Its width represents the expected statistical error 
on 20. 
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